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Abstract 

In any context in which color superconductivity arises in nature, it is 
Ukely to involve pairing between species of quarks with differing chemical 
potentials. For suitable values of the differences between chemical potentials. 
Cooper pairs with nonzero total momentum are favored, as was first realized 
by Larkin, Ovchinnikov, Pulde and Ferrell (LOFF). Condensates of this 
sort spontaneously break translational and rotational invariance, leading to 
gaps which vary periodically in a crystalline pattern. Unlike the original 
LOFF state, these crystalline quark matter condensates include both spin 
zero and spin one Cooper pairs. We explore the range of parameters for 
which crystalline color superconductivity arises in the QCD phase diagram. 
If in some shell within the quark matter core of a neutron star (or within 
a strange quark star) the quark number densities are such that crystalline 
color superconductivity arises, rotational vortices may be pinned in this shell, 
making it a locus for glitch phenomena. 



1 Overview 



The attraction between two quarks which are anti-symmetric in color renders cold 
dense quark matter unstable to the formation of quark Cooper pairs in a color 
superconducting state [0, 0, 0, |, ||, |^. If two (or more) different quark flavors 
are involved, and their Fermi momenta are the same, they pair as in the standard 
BCS state. The pairing is guaranteed because in the absence of an interaction each 
pair costs no free energy — each quark can be created at its Fermi surface — and the 
interaction then makes the system unstable against formation of a condensate of 
pairs. 

In this paper we study the situation, generic in the real world, where the Fermi 
momenta of the two species are different. If the Fermi momenta are far apart, 
no pairing between the species is possible. The transition between the BCS and 
unpaired states as the splitting between Fermi momenta increases has been studied 
in electron [|| and QCD 0, |10], |ll| superconductors, assuming that no other state 
intervenes. However, there is good reason to think that another state can occur. This 
is the "LOFF" state, first explored by Larkin and Ovchinnikov |12] and Fulde and 
Ferrell |T^ in the context of electron superconductivity in the presence of magnetic 
impurities. They found that near the unpairing transition, it is favorable to form a 
crystalline state in which the Cooper pairs have nonzero momentum. This is favored 
because it gives rise to a region of phase space where each of the two quarks in a 
pair can be close to its Fermi surface, and such pairs can be created at low cost in 
free energy. 

We study the pairing between two species whose chemical potentials differ by 
25yU and find that for a large class of interactions there is a window of 5^ within 
which states of the LOFF type are preferred over the BCS and unpaired states. 
This has important ramifications for compact star phenomenology, since it means 
that there may be a layer of crystalline quark matter inside the star. This could pin 
rotational vortices, and lead to the kind of glitch phenomena that have up to now 
been thought of as uniquely associated with the nuclear crust of neutron stars. 

In Section |^, we give a more detailed introduction to the BCS and LOFF color 
superconducting states, and their possible astrophysical applications. In Section |^, 
we describe the LOFF state in quark matter with ^ 0. We note in particular 
that, unlike in the original LOFF context, there is pairing both in J = and 
J = 1 channels. In Section ^, we derive the gap equation for the LOFF state for a 
model Hamiltonian in which the full QCD interaction is replaced by a four-fermion 
interaction with the quantum numbers of single gluon exchange. In Section |^, we 
use the gap equation to evaluate the range of 5^ within which the LOFF state 
arises. We will see that at low 5/i the translationally invariant BCS state, with 
gap Ao, is favored. At there is a first order transition to the LOFF paired 
state, which breaks translational symmetry. At pairing disappears, and 

translational symmetry is restored at a phase transition which is second order in 
mean field theory. In the weak-coupling limit, in which Aq ^ /i, we find values 
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of and Sfi2 which are in quantitative agreement with those obtained by LOFF. 
This agreement occurs only because we have chosen an interaction which is neither 
attractive nor repulsive in the J = 1 channel, making the J = 1 component of 
our LOFF condensate irrelevant in the gap equation. In Section we consider 
a more general Hamiltonian in which the couplings corresponding to electric and 
magnetic gluon exchange can be separately tuned. This leads to interactions in 
both J = and J = 1 channels, and we show how it affects the range of 6fi 
within which the LOFF state arises. In Section ^ we outline future work which 
follows immediately from what we have done and look farther ahead toward possible 
astrophysical consequences of crystalline color superconductivity. 

We recommend that the astrophysically inclined reader, interested primarily in 
the consequences of our results, read Sections ^ and |^, skipping those in between. 



2 Introduction 



2.1 Astrophysical applications of color superconductivity 



Our current understanding of the color superconducting state of quark matter leads 
us to believe that it may occur naturally in compact stars. The critical temperature 
below which quark matter is a color superconductor is generally estimated to be 
of order 10 — 50 MeV, which suggests that any quark matter which occurs within 
neutron stars that are more than a few seconds old is in a color superconducting 
state. 

This estimate of the critical temperature comes both from models whose 
parameters are tuned to reproduce zero density physics p, |, |], ||, |1^, [15|, |1^, [1^ and 
also from weak coupling methods which are quantitatively valid at asymptotically 
high densities pLIll, pOip, 
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^ , with chemical potentials 
fi ^ 10^ MeV [^. Neither class of methods can be trusted quantitatively for quark 
number chemical potentials fi ~ 400 MeV, as appropriate for the quark matter 
which may occur in the cores of neutron stars. Still, both methods agree that the 
gaps at the Fermi surface are of order tens to 100 MeV, with critical temperatures 
about half as large. 

It is therefore important to look for astrophysical consequences of color 
superconductivity. As a Fermi surface phenomenon, it has little effect on the 
equation of state, and hence little effect on the radius of a compact star. There are 
nevertheless several effects of color superconductivity under active investigation. 
The color superconductivity of quark matter in neutron stars influences the 
evolution of magnetic fields within the quark matter (see also ||31[]). Cooling 
by neutrino emission is also affected (see also [^). In quark stars, the 
physics of the instability to r-mode oscillations is dramatically affected by color 
superconductivity although this is not the case for neutron stars with quark 
matter present only in their cores |3^, |3^. Furthermore, the phase transition at 
which color superconductivity sets in as a hot proto-neutron star cools may yield 
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a detectable signature in the neutrinos received from a supernova Finally, 
one goal of the present paper is to motivate an investigation of the possibility that 
(some) pulsar glitches may originate in quark matter. 

If two species of fermion experience an attractive interaction, and their Fermi 
momenta are the same, they pair in the standard BCS state. The pairing is 
guaranteed because in the absence of an interaction each pair costs no free energy 
(each quark can be created at its Fermi surface) , and the interaction then makes the 
system unstable against formation of a condensate of pairs. In the QCD context, if 
there are two flavors of quarks with equal Fermi momenta, quarks of two colors and 
two flavors pair p, ^ while if there are three flavors of quarks, all nine quarks pair in 
a pattern which locks color and flavor symmetries, breaking chiral symmetry |3^ . 
These idealizations are very instructive, but in any physical context, the up, down 
and strange quarks will all have different Fermi momenta. To give the reader some 



sense for typical scales in the problem, we give an illustrative example |^|. The 



s 



numbers in this paragraph assume that the quarks are noninteracting fermions — 
clearly a bad assumption — and so should certainly not be construed as precise. 
Consider quark matter with average quark chemical potential /i = 400 MeV, made 
of massless up and down quarks and strange quarks with mass = 300 MeV. (M, 
is a density dependent effective mass; this adds to the uncertainty in its value.) If 
the strange quark were massless, quark matter consisting of equal parts u, d and 
s would be electrically neutral. In our illustrative example, on the other hand, 
electric neutrality requires a nonzero density of electrons, with chemical potential 
/^e = 53 MeV. Charge neutrality combined with the requirement that the weak 
interactions are in equilibrium determine all the chemical potentials and Fermi 
momenta: 

/i« = - |/Ue = 365 MeV, pp = 

/Xd = /i + |/ie = 418 MeV, = fid, 

^ (2.1) 

/X, = /i + i/ie = 418 MeV, p'p = ^/x2 _ ^2 = 290 MeV, 



fie = 53 MeV, PF = f^ 



e 



The baryon number density pb = (l/37r^) [(p^)^ + {ppY + (pp)^] is 4 times nuclear 
matter density.^ As one goes deeper into a neutron star, /i increases, Mg decreases 
somewhat, and fi^ and all differences between the quark Fermi momenta decrease.^ 
In this paper, we investigate the consequences of pairing between quarks with 

^Had we chosen = 200 MeV, we would have obtained fi^ = 24 MeV, p}. = 384 MeV, 
Pp = 408 MeV, pp = 356 MeV and a, pB of 5 times nuclear matter density. 

^Note that in a neutron star with a quark matter core, regions of purely hadronic and purely 
quark matter are separated by a mixed phase in which neither the hadronic regions nor the quark 
matter regions are separately charge neutral |Q. The electrically charged quark matter in these 
regions will have Fermi momenta which differ qualitatively from those in our example, with pp 
less than either pp or pp . 
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differing Fermi momenta. For simplicity, we restrict our explicit calculations to 

the case of two massless quarks with differing chemical potentials and fi^, which 
we write as 

IJ,d = ft + 5n; /i„ = /i-(5/i. (2.2) 

We expect similar phenomena to those we describe to arise wherever any one of 
— Pf\ or — Pf\ \Pf ~ Pf\ falls within a suitable range, but we leave the 
investigation of quark matter with u, d, and massive s quarks to future work. We 
also work at zero temperature throughout. 

2.2 Isotropic (non-LOFF) pairing 

In the color superconducting phase for two massless quark flavors at the same 
chemical potential /i, the condensate consists of quark-quark pairs which are 
flavor singlets and color 3 antitriplets (and hence also spin singlets, to obey Pauli 
statistics). Pairing is of the BCS type: a red up quark of momentum p pairs 
with a green down quark of momentum — p of the same helicity, so that the spins 
are antiparallel. The blue quarks are left unpaired. Such pairing is strongest in the 
vicinity of the Fermi surface, for | |p| — /i| < Aq, where Aq is the BCS gap parameter. 

If, instead, the Fermi momenta are sufficiently different, no BCS pairing is 
possible. It is no longer possible to guarantee that the formation of pairs lowers 
the free energy, because in the BCS state the two fermions in a pair have equal 
and opposite momentum, so at most one member of each pair can be created at its 
Fermi surface. The other member costs non-zero free energy, which the attractive 
interaction may be unable to compensate. 

Assuming that no other state intervenes between the BCS state and the state 
with no condensate, we can apply the results first derived by Clogston and 
Chandrasekhar in the context of pairing between spin-up and spin-down electrons 
with differing Fermi momenta ||^. For small enough the favored BCS state has 
coincident Fermi surfaces, p% = p% = Ji because this maximizes the pairing and thus 
the gain in interaction energy.^ We denote the gap in this BCS state by Aq. The 
free energy of this BCS state must be compared to that of the unpaired or "normal" 
state in which the quarks simply distribute themselves in Fermi seas with p^ = fiu, 
Pf = A*rf and no condensate forms. The BCS state is the stable ground state of 
the system only when its negative interaction energy offsets the large positive free 
energy cost associated with forcing the Fermi seas to deviate from their normal state 
distributions. If Aq ^ /i and 5/i <^ ft, the free energy of the BCS state relative to 
that of the normal state at a given 6fi is 

i^BCS - i^normal = ^ - . (2.3) 

■^If one tries to construct a "BCS-like" state which has = /i„ and pf, = fXd and consequently 
no pairing for < p < p%, one finds ||3^ that this state has a higher free energy than the BCS 
state (in which p|n = pf = p.). The gain in free energy associated with choosing pp = ^„ and 
Pf = does not compensate for the fost pairing energy. 
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The coefficient I/Stt^ depends on the number of fermion species which pair and is 
appropriate to the case of interest to us. Clogston and Chandrasekhar concluded 
that the BCS state is favored for 5/i < = Ao/\/2. (The relation = Aq/\/2 
is exact only in the weak-coupling limit in which Aq <^ fi.) At 6fj, = Sfii, there is a 
first order phase transition at which the gap parameter drops discontinuously from 
Aq to zero: for 6fi < 6fii, the system is in the BCS phase, unperturbed from the 
the 6fi = state. 

This analysis is modified in an interesting way at nonzero temperature, as was 
discussed by Lombardo and Sedrakian in the context of pairing between neutrons 



and protons in nuclei |^0[. Thermal excitations smear out the normal state Fermi 
surfaces, making pairing between thermally excited states above the lower Fermi 
surface and below the upper Fermi surface possible. As a consequence, as T is 
increased from zero, there is a range of T within which 6fii is larger than at T = 0. 
At still higher temperatures, of course, all pairing is lost. 

In applying the work of Clogston and Chandrasekhar to color superconductivity, 
there have been two extensions to their analysis. First, recall that only two colors 
of up and down quarks pair. In describing the mixed phase associated with the first 
order phase transition, one must take careful account of the unpaired blue quarks. 



This has been done by Bedaque ]TT[. Second, in Refs. |T0| the transition between 
the color-flavor locked phase and the two-flavor color superconducting phase has 
been studied, under the assumption that yUn = A^d = /^s but with pp ^ p^'^ because 
of the nonzero strange quark mass. The ffist order transition that these authors 
describe is similar to that of Clogston and Chandrasekhar, as it is associated with 
the unpairing of us and ud Cooper pairs, but it differs in that the analogue of the 
normal state is one in which u and d quarks remain paired. As we have seen above, 
treating a realistic situation requires relaxing the assumption of equal chemical 
potentials. 

2.3 Non-isotropic (LOFF) pairing 

The Clogston and Chandrasekhar analysis of the ffist order unpairing transition 
assumes that the only possible phases of the system are a BCS phase and the 
normal phase. However, there is good reason to think that another state can occur 
in the crossover region between BCS and no pairing. As was first realized by Larkin 



and Ovchinnikov |T2[ and Fulde and Ferrell |T3| (LOFF), whereas the BCS state 
requires pairing between fermions with equal and opposite momenta, when 6fi ~ 6fii 
it may be more favorable to form a condensate of Cooper pairs with nonzero total 
momentum. By pairing quarks with momenta which are not equal and opposite, 
some Cooper pairs are allowed to have both the up and the down quarks on their 
respective Fermi surfaces even when 6fi ^ 0. LOFF found that within a range of Sn 
near Sfi^, a condensate of Cooper pairs with momenta q + p and q — p (see Figure |ID 
is favored over either the BCS condensate or the normal state. Here, our notation is 
such that p specifies a particular Cooper pair, while q is a fixed vector, the same for 
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all pairs, which characterizes a given LOFF state. The magnitude |q| is determined 
by minimizing the free energy; the direction of q is chosen spontaneously. The 
resulting LOFF state breaks translational and rotational invariance. In position 
space, it describes a condensate which varies as a plane wave with wave vector 2q. 

Once one has demonstrated an instability to the formation of a plane wave, 
it is natural to expect that the state which actually develops has a crystalline 
structure. Larkin and Ovchinnikov in fact argue that the favored configuration is a 
crystalline condensate which varies in space like a one-dimensional standing wave, 
cos(2q-r). Such a condensate vanishes along nodal planes |jl2|. Subsequent analyses 
suggest that the crystal structure may be more complicated. Shimahara has 
shown that in two dimensions, the LOFF state favors different crystal structures at 
different temperatures: a hexagonal crystal at low temperatures, square at higher 
temperatures, then a triangular crystal and finally a one- dimensional standing wave 
as Larkin and Ovchinnikov suggested at temperatures that are higher still. In three 



dimensions, the question of which crystal structure is favored seems unresolved |42 



LOFF did their analysis in the same context as that of Clogston and 
Chandrasekhar: electromagnetic superconductivity in a magnetic field which causes 
a Zeeman splitting while not inducing screening currents. They were seeking to 
model the physics of magnetic impurities in a superconductor. Magnetic effects on 



the motion of the electrons and the scattering of electrons off non-magnetic 
impurities ^ disfavor the LOFF state. Although signs of the BCS to LOFF 
transition in the heavy fermion superconductor UPd2Al3 have been reported 



the interpretation of these experiments is not unambiguous |^7[. It has also been 
suggested that the LOFF phase may be more easily realized in condensed matter 



systems which are two-dimensional |^ , |48[| or one- dimensional |^ , both because in 



these cases 5/i2 is larger than in three-dimensional systems and because the magnetic 
field applied precisely parallel to a one- or two-dimensional system does not affect the 
motion of electrons therein. Evidence for a LOFF phase in a quasi-two-dimensional 
layered organic superconductor has recently been reported |^ . 

None of the difficulties which have beset attempts to realize the LOFF phase 
in a system of electrons in a magnetic field arise in the QCD context of interest 
to us. Differences between quark chemical potentials are generic and the physics 
which leads to these differences has nothing to do with the motion of the quarks. 
We therefore expect the original analysis of LOFF (without the later complications 
added in order to treat the difficulties in the condensed matter physics context) 
to be a good starting point. In this paper we use an analysis based on that 



originally done by Fulde and Ferrell |]T3|, but described in more detail by Takada 



and Izuyama [^|, to argue that for appropriate values of quark number densities, 
the color superconducting gap may vary periodically in space, forming a crystalline 
pattern. More precisely, what we will demonstrate is that if some difference between 
chemical potentials falls in the appropriate range, quark matter is unstable to the 
spontaneous breaking of translational invariance by the formation of condensates 
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Figure 1: The momenta ku and k,^ of the two members of a LOFF-state Cooper pair. 
We choose the vector q, common to all Cooper pairs, to coincide with the 2;-axis. 
The angles and ad{p) indicate the polar angles of k„ and k^, respectively. 



which vary in space like a plane wave. Following Larkin and Ovchinnikov ||T2 
we expect that once there is an instability to the formation of plane waves the 
condensate that results will be crystalline, but we leave the determination of the 
crystal structure of the condensate to future work. 



3 The LOFF state 



We begin our analysis of a LOFF state for quark matter by constructing a variational 
ansatz for the LOFF wavefunction. As motivated by the preceding discussion, we 
consider Cooper pairs which consist of an up quark and a down quark with respective 
momenta 

k« = q + P, kd = q-p, (3.1) 

so that p identifies a particular quark pair, and every quark pair in the condensate 
has the same nonzero total momentum 2q. This arrangement is shown in Figure |l|. 
The helicity and color structure are obtained by analogy with the "2SC" state as 
described in previous work p, Q : the quark pairs will be color 3 antitriplets, and in 
our ansatz we consider only pairing between quarks of the same helicity. 

With this in mind, here is a suitable trial wavefunction for the LOFF state with 



7 



wavevector q [0 |T§ 0: 
|vl/q) = 44|0), 

4 = n (cos^l(p) + 6"^V«-(P) sin^i(p) aL.(q+p) aU(q-p)) 

^^^'"•^ (3.2) 

X n 4«a(q+p) X n 4d/3(q-p)' 

= as above, L ^ R, 

where a, j3 are color indices, m, and L, i? are the usual flavor and helicity labels, 
and is the particle creation operator (for example, a^^^ creates a left-handed 
down quark with color a). The ^'s and ^'s are the variational parameters of our 
ansatz: they are to be chosen to minimize the free energy of the LOFF state, as 
described in the next section. The first product in equation ( p.2| ) creates quark 
pairs within a restricted region V of the total phase space. This allowed "pairing 
region" will be discussed below. The next product fills a "blocking region" Eu with 
unpaired up quarks: these are up quarks with momenta q + p for which there are 
no corresponding down quarks with momenta q — p. The final product fills the 
blocking region Ed with unpaired down quarks. The ansatz does not contain a term 
that would create antiparticle pairs: we have checked the effect of such a term and 
found that it has no qualitative effect on our results. 

To complete the specification of our ansatz we need to describe the allowed 
pairing and blocking regions in phase space. These regions are largely determined 
by Pauli blocking as a result of populated Fermi seas. In the absence of pairing 
interactions, the system is in the "normal" state and up and down quarks are 
distributed in Fermi seas with Fermi momenta = /i„ and pf, = /x^, respectively 
(recall that we consider massless quarks only, so the single particle energy of a quark 
with momentum k is e(k) = |k|). An up quark carries momentum k^j = p + q; 
in p-space, therefore, the Fermi sea of up quarks corresponds to a sphere of radius 
Aiji = fi—Sfi centered at — q. Similarly, a down quark carries momentum = — p+q, 
giving a sphere in p-space of radius /x^ = /i + 5/i centered at +q. The two offset 
spheres are shown in Figure ^ (we have drawn the case |q| > 6^ so that the two 
Fermi surfaces intersect in p-space). In the limit of arbitrarily weak interactions, the 
blocking region Eu corresponds to the lower shaded area in the figure: pairing does 
not occur here since the region is inside the Fermi sea of up quarks, but outside the 
Fermi sea of down quarks. Similarly the upper shaded area is the blocking region Ed- 
The entire unshaded area is the pairing region V: it includes the region inside both 
spheres, where hole-hole pairing can occur, and the region outside both spheres, 
where particle-particle pairing can occur. 

We can now explain how the LOFF wavefunction ansatz can describe the normal 
state with no condensate: we choose Ol{p) = ^_r(p) = vr/2 for p inside both Fermi 
spheres, and otherwise all the ^'s are zero. With this choice the first term in 
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Figure 2: The LOFF phase space, as a function of p (Eq. ( |3.1| )). We show the 
Py = plane, (a) The phase space in the hmit of arbitrarily weak interactions. In 
the shaded blocking regions and Bd, no pairing is possible. In the inner unshaded 
region, an interaction can induce hole-hole pairs. In the outer unshaded region, an 
interaction can induce particle-particle pairs. The region V (Eq. ( ^.2D ) is the whole 
unshaded area, (b) When the effects of interactions and the formation of the LOFF 
state are taken into account, the blocking regions shrink. The BCS singularity 
occurs on the dashed ellipse, defined by e„ + = -|- /x^, where making a Cooper 
pair costs no free energy in the free case. 

Eq. (|3.2| ) fills that part of each Fermi sea corresponding to the inner unshaded 
region of Figure The B^ and Bd terms fill out the remainder of each Fermi sea 
to obtain the normal state. Note that in the absence of pairing, the normal state 
can be described with any choice of q. The most convenient choice is q = 0, in 
which case k„ = = p, Bu vanishes, and Bd is a spherical shell. Other choices 
of q correspond to choosing different origins of k„-space and k^-space, but in the 
absence of any interactions this has no consequence. Once we turn on interactions 
and allow pairing, we expect a particular |q| to be favored. 

The phase space picture changes slightly when pairing interactions are included: 
the blocking regions are smaller when a LOFF condensate is present, as indicated 
in Figure ^d. We will account for this effect in the next section. With smaller 
blocking regions, a larger portion of the phase space becomes available for LOFF 
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pairing. Such pairing is guaranteed to be energetically favorable when it costs zero 
free energy to create an up quark and a down quark, since these quarks can then 
pair to obtain a negative interaction energy. The zero free energy condition is 



where e(k) is the single particle energy of a quark with momentum k. For massless 
quarks, we obtain |q + p| + |q — p| = 2/i, which describes an ellipsoidal surface in 
p-space. This surface is indicated by the ellipse shown in Figure ^b; notice that the 
ellipsoid and the two Fermi surfaces all intersect at a circle. 

If the interaction is weak, we expect LOFF pairing to be favored in a thin layer 
of phase space around this ellipsoid. This is manifest in the gap equation derived in 
the next section (Eq. (|4.9|) ) in which, as in BCS theory, we find a divergent integrand 
on this ellipsoid in the absence of pairing. Pairing smoothes the divergence. As the 
interaction gets stronger, the layer of favored pairing gets thicker. If there were 
no blocking regions, we could use the entire ellipsoid, just as BCS pairs condense 
over the entire spherical surface |p| = fi in the symmetric, S/i = |q| = case. 
However, as shown in Figure ^3, the blocking regions exclude pairing over most of 
the ellipsoid, leaving a ribbon of unsuppressed LOFF pairing in the vicinity of the 
circle where the Fermi surfaces intersect. This agrees with our expectation for the 
particle distribution in the LOFF state: it is as in the normal state, except that 
there is a restricted region (around the aforementioned ribbon) where each quark in 
a pair can be near its Fermi surface. 

Although the constant single-particle energy contours for noninteracting up and 
down quarks cross in p-space (see Figure ^), we emphasize that the Fermi surfaces 
of up and down quarks do not cross in momentum (k„- and k^-) space. The p-space 
ribbon of unsuppressed pairing corresponds to unsuppressed pairing between up and 
down quarks with momenta around k-space ribbons near their respective (disjoint) 
Fermi surfaces. 

In the limit of arbitrarily weak interactions, the ribbon in momentum space along 
which pairing is unsuppressed shrinks, as the blocking regions grow to exclude all of 
the ellipsoid except the one-dimensional circle at which the two spheres in Figure ^ 
intersect. This circle has insufficient phase space to lead to a singularity in the gap 
equation: the integrand is singular on this circle, but the integral does not diverge. 
Therefore, the LOFF state is not guaranteed to occur if one takes the weak coupling 
limit at fixed 5fi. In this respect, the LOFF state is like the BCS state at nonzero 
6fi: for weak coupling, Aq and because the BCS state can only exist if it has 
Ao > v^(5/x, it must vanish for couplings weaker than some threshold. We shall see, 
however, that at any fixed weak coupling, the LOFF state, like the BCS state, is 
guaranteed to occur at some 5/i: the BCS state arises if 6fi < 6fii and the LOFF 
state arises if 6fii < 6fi < 5^2- 

One of the most striking features of the LOFF state is the spin structure of 
the condensate. The familiar "2SC" state pairs quarks of the same helicity and 




(3.3) 
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opposite momentum, so the spins are antiparallel and the quarks are arranged in 
an antisymmetric combination to form spin singlet Cooper pairs. The LOFF state 
also pairs quarks of the same helicity, but now the quark momenta are no longer 
antiparallel, as can be seen from Figure |l|. Therefore the LOFF Cooper pairs cannot 
be spin singlets: they are superpositions of both spin zero and spin one. This is 
revealed explicitly by evaluating the nonzero {ipip) expectation values in the LOFF 
state: 

where j are flavor indices (1 = up, 2 = down), a, (3 are color indices, C = ^7^7^, 
L = (1 — 75)72 is the usual left-handed projection operator, and a^u = (i/2)[7^, 7,^]. 
The constants and Fg are left-handed J = and J = 1 condensates, respectively. 
F^ and are defined analogously. The F's can be expressed in terms of the 
variational parameters of the LOFF wavefunction: 



F^ = — sm6'i(p) cos0i(p)e"^^^P^ sm(^- 



g-*0{p) 



pG-P 

4 



= ^j:sn.eMcoseMe^^^^-^ sin(^MPA^iMPl)e- 



■i(/>(p) 



(3.5) 



Here V is the spatial volume of the system, a(p) are the polar angles of the quark 
momenta, as in Figure |l], and the dependence on the azimuthal angle (p follows from 
our use of the spinor conventions described in Refs. 0, |^, The expressions for 
F^ and are the same as those in (|3.5| ) except that <f){p) is replaced by vr — 0(p). 
In Eq. ( |3.5| ) and throughout, (1/V") becomes / cPp/{27iY in an infinite system. 

Once we have derived a gap equation by minimizing the free energy with respect 
to these variational parameters, we expect the condensates to be simply related to 
gap parameters occurring in the gap equation. We will see explicitly how F^ and 
F^ are determined in the next section. 

Notice that the condensates of Eq. ( p.4|) are plane waves in position space by 
virtue of the nonzero momentum 2q of a Cooper pair. describes pairing which 
is antisymmetric in color, spin, and flavor, while F^ describes pairing which is 
antisymmetric in color but symmetric in spin and flavor (in each case, Pauli statistics 
are obeyed). In the original LOFF condensate of electrons there can be no F^, since 
electrons have no color or flavor, so that only the spin antisymmetric pairing is 
possible. 

The J = condensates (ipCLip), (ipCRip) are Lorentz scalars (mixed under 
parity), while the J = 1 condensates (ipCLa^^ip) y {'ipCRa^^ip) are 3- vectors (also 
mixed under parity) which point in the z-direction, parallel to the spontaneously 
chosen direction q of the LOFF state. Because the ansatz contains a J = 1 
component, it would be interesting to generalize it to include the possibility of 
LR pairing, in addition to LL and RR pairing. We discuss this further in Section 
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The possibility of a LOFF phase in QCD has been mentioned briefly in a different 
context. In their analysis of quark matter with a very large isospin density (with 
large Fermi momenta for down and anti-np quarks) Son and Stephanov have noted 
that if the d and u Fermi momenta differ suitably, a LOFF phase will arise . 



In the physically realizable context of large baryon number density, pairing 



between quarks and holes with nonzero total momentum has also been discussed |53 



5H ^ p^ . This results in a condensate with the quantum numbers of (qq), which 
varies in space with a wave number equal to in contrast, the LOFF phase 
describes a diquark condensate which varies with a wave number 2|q| comparable 
to 2(5/i. The crystalline chiral condensate ||5^ is favored in QCD at asymptotically 
high densities only if the number of colors is very large [^|, greater than about 



Nc = 1000 [0, 0. At lower densities, where the interaction is stronger, the 
crystalline chiral condensate may arise in QCD with fewer colors [^. Apparently, 
however, in QCD with Nc = 3 this phase is not favored (although it is close to 
being competitive) even when the coupling is so large that Ao//i > 1/2. Note 
that crystalline color superconductivity is guaranteed to occur at arbitrarily weak 
coupling for suitably chosen 6fi, while a crystalline chiral condensate cannot form 
anywhere in the phase diagram if the coupling is weak. 



4 The gap equation and free energy 

Having presented a trial wavefunction for the LOFF state, we now proceed to 
minimize the expectation value of the free energy (F) with respect to the variational 
parameters of the wavefunction (the ^'s and ^'s of equation ( p.2|) ) to obtain a LOFF 
gap equation. The free energy is F = H — fiuN^— fidNd, where H is the Hamiltonian, 
and Nu and A^^^ are the number operators for up and down quarks, respectively. We 
choose a model Hamiltonian which has a free quark term Hq and an interaction term 
Hj, and write the free energy a.s F = Fo + Hj, where Fq = Hq — fiuNu — fidNd is the 
free energy for noninteracting quarks. To describe the pairing interaction between 
quarks, we use an NJL model consisting of a four-fermion interaction with the color 
and flavor structure of one-gluon exchange: 

Hi = lj d'x [G^(^7°T^^)(^7°T^^) - GM(^yT^^)(^yT^^)] (4-1) 

where the are the color SU{3) generators, normalized so that tr(T"^T'^) = 25^^. 
Notice that we have relaxed some constraints on the spin structure of one-gluon 
exchange: we allow for the possibility of independent couplings Ge and Gm for 
electric and magnetic gluons, respectively. This spoils Lorentz boost invariance but 
there is no reason to insist on boost invariance in a finite-density system. Indeed, in 
high density quark matter we expect screening of electric gluons but only Landau 
damping of magnetic gluons, and we might choose to model these effects by setting 
Ge ^ Gm- We postpone a discussion of these issues and their implications for the 
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LOFF state until Section |^. For now, we restrict ourselves to the case of Lorentz 
invariant single gluon exchange, by letting Ge = Gm = G > 0. 

We need to evaluate the expectation value of F in the LOFF state to obtain an 
expression for the free energy of the system in terms of the variational parameters 
of the ansatz. The noninteracting part of the free energy is simply 

(Fo) = 5^ 2(|q + p| - + 5^ 2(|q - p| - /i,) 
peBu p&Bd 

+ ^ 2(|q + p| + |q - p| - /i„ - fid) sin^ ^^(p) (4.2) 
pev 

+ (same, with L ^ R). 

The first and second terms represent the contributions of the unpaired left-handed 
up and down quarks, respectively. The third term gives the (noninteracting) free 
energy of the left-handed quark pairs. The three terms are all repeated with L 
replaced by R to include the free energy for the right-handed quarks. The factors of 
two in equation ( [4.2| ) appear because there are two quark colors ( "red" and "green" ) 
involved in the the condensate. The "blue" quarks do not participate in the pairing 
interaction and instead behave as free particles: the blue up and down quarks fill 
Fermi seas with Fermi momenta = /i^ and p% = fid, respectively. Below, we 
will want to compare the free energy of the LOFF, BCS and normal states. Since 
at any given /i„ and fid the free energy of the spectator quarks is the same in all 
three states, we can neglect these blue quarks in the remainder of our analysis even 
though they do contribute to the total free energy. 

The expectation value of Hj gives the total binding energy of the pairing 
interaction: 

m = -IGV {\T^^\' + (4.3) 

where the F^'s are the J = LOFF condensates defined in equations (|3.5|) . These 
condensates are simply related to J = LOFF gap parameters defined as 

The gap parameters correspond to IPI Green's functions and are the quantities 
which will appear in the quasiparticle dispersion relations and for which we will 
derive the self-consistency conditions conventionally called gap equations. We see 
from Eq. ( |4.3| ) that with G > the interaction is attractive in the J = channel 
and is neither attractive nor repulsive in the J = 1 channel. 

Our ansatz breaks rotational invariance, so once J = pairing occurs (Ta 7^ 0) 
we expect that there will also be J = 1 pairing (F^ 7^ 0). As we have seen, this 
arises even in the absence of any interaction in the J = 1 channel as a consequence 
of the fact that the momenta of two quarks in a Cooper pair are not anti-parallel 
if q 7^ 0. Because (H) is independent of F^, the quasiparticle dispersion relations 
must also be independent of F^. That is, the J = 1 gap parameter must vanish: 
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= 0. In Section we shall see by direct calculation that A^ is proportional 
to {Ge — Gm)^b- In the present analysis with Ge = Gm, therefore, A^ = while 

The ^'s are chosen to cancel the azimuthal phases 0(p) in equations ( p.5|) . By 
this choice we obtain maximum coherence in the sums over p, giving the largest 
possible magnitudes for the condensates and gap parameters. We have 

^l(p) = 0(p) + ¥'L, ^^(P) = - 0(p) + </?ij (4.5) 

where (pi and (pn are arbitrary p-independent angles. These constant phases do not 
affect the free energy — they correspond to the Goldstone bosons for the broken 
left-handed and right-handed baryon number symmetries — and are therefore not 
constrained by the variational procedure. For convenience, we set cpL = <Pr = and 
obtain condensates and gap parameters that are purely real. 

The relative phase (fL — fR determines how the LOFF condensate transforms 
under a parity transformation. Its value determines whether the J = condensate 
is scalar, pseudoscalar, or an arbitrary combination of the two and whether the 
J = 1 condensate is vector, pseudovector, or an arbitrary combination. Because 
single gluon exchange cannot change the handedness of a massless quark, the left- 
and right-handed condensates in the LOFF phase are not coupled in the free energy 
of Eq. (|4.3| .) Our choice of interaction Hamiltonian therefore allows an arbitrary 
choice of ipi — ^r- A global U{1)a transformation changes </?l — ^r, and indeed 
this is a symmetry of our toy model. If we included ?7(l)A-breaking interactions 
in our Hamiltonian, to obtain a more complete description of QCD, we would find 
that the free energy depends on ip^ — ipR, and thus selects a preferred value. For 
example, had we taken Hj to be the two-flavor instanton interaction as in Ref. 0, 
the interaction energy would appear as r'^*r^ -|- r-^F-^* instead of as in ( |4.3| ). This 
would enforce a fixed phase relation (f^ — fR = 0, favoring condensates which are 
parity conserving P, ^ . 

We now apply the variational method to determine the angles ^(p) in our trial 
wavefunction, by requiring that the free energy is minimized: d{F) /d6{p) = 0. This 
is complicated by the fact that the pairing region V and the blocking regions Bu 
and Bd are themselves implicitly dependent on the 6 angles: these angles determine 
the extent of the LOFF pairing, and the phase space regions V, Bu and Bd change 
when a condensate is present, as mentioned in Section ^ For now we simply ignore 
any 6'-dependence of the phase space regions; our result will nevertheless turn out 
to be correct. Everything is the same for left and right condensates so we hereafter 
drop the L and R labels. Upon variation with respect to 6'(p), we obtain 

ta„2«(p) = ^^^«^l/H^ (4.6) 

|q + p| + |q - p| - /i« - yUd 

where /9yi(p) = a?i(p) + ^^(p) is the angle between the two quark momenta in a 
LOFF pair, as shown in Figure 0. Notice that the denominator on the right hand 
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side of the above expression vanishes along the elhpsoidal surface of optimal LOFF 
pairing described in Section |^. When q = 0, the quark momenta are antiparallel so 
/3yi(p) = TT and Eq. ( |4.6| ) reduces to the simple BCS result: tan26' = Ayi/(|p| — fl). 

With the 6 angles now expressed in terms of a gap parameter A^, we turn to 
the LOFF quasiparticle dispersion relations. They can be obtained by taking the 
absolute value of the expressions 



Ei{p) = 5/i + i(|q + p| - |q-p|) 



+ ly/{h + P\ + \^-p\- 2/^)2 + AAl sin2(i/3^(p)) 

^2(p) = - i(|q + p| - |q - p|) 

+ + Pl + lq - Pl - + 4Ai sin2(i/3^(p)) 



(4.7) 



whose meaning we now describe. For regions of p-space which are well outside both 
Fermi surfaces, Ei {E2) is the free energy cost of removing a LOFF pair and adding 
an up quark with momentum q + p (a down quark with momentum q — p) . For 
regions of p-space which are well inside both Fermi surfaces, Ei {E2) is the free 
energy cost of removing a LOFF hole pair and adding a down hole with momentum 
q — p (an up hole with momentum q + p). Where the Fermi surfaces cross in p- 
space and pairing is maximal, both quasiparticles are equal superpositions of up and 
down. In the region of p-space which is well inside the up Fermi surface but well 
outside the down Fermi surface, Ei is negative, corresponding to a domain in which 
it is energetically favorable to have an unpaired up quark with momentum q + p 
rather than a (q+p, q — p) quark pair. Similarly, E2 is negative where it is favorable 
to have an unpaired down quark with momentum q — p rather than a LOFF pair. 
Equations ( [4.7|) allow us to finally complete our description of the LOFF phase by 
specifying the definitions of the phase space regions V, Bu and B^- The blocking 
region Bu is the region where Ei{p) is negative, and unpaired up quarks are favored 
over LOFF pairs. Similarly Bd is the region where E2{p) is negative. The regions 
ii^i < and E2 < are shown as the shaded areas in Figure ^ for A a = 0, and 
in Figure 0b for 7^ 0. LOFF pairing occurs in the region where Ei and E2 are 
both positive: 

P = {p|Ei(p) > and E2(p) > 0} (4.8) 

corresponding to the entire unshaded regions of Figure 0. The actual quasiparticle 
dispersion functions are |ii^i(p)| and |i?2(p)h they are nonnegative everywhere, since 
they represent energies of perturbations of the LOFF state which is the presumed 
ground state of the system.0 In the blocking regions, elementary excitations are 
created by replacing an unpaired quark with a quark pair, and vice versa in the 



^ Since the LOFF condensate contains pairs with momentum 2q, the momentum of its 
quasiparticle excitations is only defined modulo 2q. The momentum, modulo 2q, of a quasiparticle 
of energy |i^i(p)| is p mod 2q. 
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pairing region. When q = 0, Eqs. (O) reduce to the more famihar BCS resuh: 



With the boundaries of the blocking regions specified, one can verify by exphcit 
calculation that the variation of these boundaries upon variation of the ^'s does not 
change the free energy. This can be understood as follows. Notice that because we 
can create zero-energy quasiparticles on the boundaries of the blocking regions, there 
is no actual energy gap in the excitation spectrum of the LOFF state. The change 
in (F) due to variation of the boundaries of the blocking regions is zero because this 
variation simply creates zero-free-energy quasiparticles on these boundaries. This 
justifies our neglect of the ^-dependence of the phase space regions in the derivation 
of Eq. (pl). 

Substituting the expression (|4.6|) for the 9 angles into the expression for 



the Ta condensate, and using the relation = GTa, we obtain a self-consistency 
equation for the gap parameter A^^: 

^^^y 2sin^(l/?^(p)) ^^g^ 

^ PGP A/(|q + Pl + |q - Pl - + 4Ai sin2(i/?^(p)) ' 



This can be compared to the BCS gap equation, obtained upon setting q = and 
eliminating the blocking regions: 

2Gv^ 1 , , 

(4.10) 



Note that in the LOFF gap equation ( |4.9| ), the gap parameter appears on the right 
hand side both explicitly in the denominator and also implicitly in the definition 
of the pairing region V, as given in ( [4.8| ). This means that if the q — limit is 
taken at fixed 6fi, the LOFF gap equation will only become the BCS gap equation 
if the blocking regions vanish in this limit. This happens if, as q ^ 0, A^ tends 
to a limiting value which is greater than 6fi. A state with A^ < 6fi and q = 
is "BCS-like", in that the Cooper pairs have zero momentum, but has no pairing 
within a region < |p| < pp. Such states always have higher free energy than the 
BCS state obtained simply by solving the gap equation ([4.10|) , appropriate if there 



are no blocking regions and p^ = pp ||39|| . 



In the next section we will solve the LOFF gap equation ( [4.9|) and determine 
the circumstances in which the LOFF state is the true ground state of the system. 
Once we have obtained a solution to the gap equation (|4.9| ) for A^, the condensates 
are given by = A^/G and 



2 sr^ 2AAsin(i/5A(p))sin(i/5B(p)) 



B 



V 



PGP x/dq + p| + |q - p| - 2/i)2 + 4Ai sin2(i/3^(p)) 



where /9_b(p) = au(p) — ad(p)- (See Figure |l|.) We now see explicitly that if the 
interaction is attractive in the J = channel, creating a nonzero F^ and A^, a 
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nonzero J = 1 condensate is induced regardless of the fact that there is no 
interaction in the J = 1 channeL As a check, note that if q = 0, sin(i/5A(p)) = 1 
and sin(i/?B(p)) is given by the cosine of the polar angle of p. The right hand side 
of ( [4.11|) therefore vanishes upon integration, and Tb vanishes when q = as it 
should. It is now apparent that two features contribute to a nonzero r^. The first 
is that the momenta in a quark pair are not antiparallel, which leads to the factors 
of sin(i/5^(p)) in Eq. (|4.11|) . The second is that the pairing region is anisotropic, 
since if it were not the factor of sin(|/3B(p)) would ensure that the right-hand side 
of ( [4 .111) vanishes upon integration. 

As written, the gap equations ( |4.9| ) and ( [4.10| ) are ultraviolet divergent. In QCD, 
of course, asymptotic freedom implies that the interaction between quarks decreases 
at large momentum transfer and we have not yet represented this fact in our toy 
model. In previous work j^, ^ ^ , we chose to mimic the effects of asymptotic freedom 
(and to render the right hand side of the gap equation finite) by introducing a form 
factor associated with each fermion leg in the four-fermion interaction. This is not 
a good strategy when q 7^ 0. The two incident quarks carry momenta q + p and 
q — p while the outgoing quarks carry momenta q + p' and q — p'. Were we to cut 
off these four momenta with form factors on each leg, we would have a cutoff which 
depends explicitly on q. This is not a good representation of what happens in full 
QCD, in which the condition for when the interaction becomes weak is determined 
by the momentum p — p' transferred through the gluon and has nothing to do with 
q. For simplicity, we choose to introduce a hard cutoff in our NJL model, rather 
than a smooth form factor, and choose simply to cut off the momentum p. This is 
not equivalent to cutting off the momentum transfer, but has the desired feature of 
being a q-independent cutoff. That is, we limit the integration region to |p| < A in 
the BCS gap equation ( [4.1(]| ) and to {p G P and |p| < A} in the LOFF gap equation 
(|4.9| ). In the BCS case, this criterion is equivalent to cutting off the momentum of 
each fermion leg. In the LOFF case, it is not equivalent and is more appropriate. 
The choice we have made is not the only q-independent cutoff one might try. For 
example, we have also obtained results upon cutting off momenta outside a large 
ellipsoid in p-space, confocal with the centers of the two Fermi spheres in Figure |^, 
but have found that this makes little difference relative to the simpler choice of the 
large sphere |p| < A. 



5 Results 

We solve the gap equation ([4.9|) numerically (and analytically in the limit ^ 
6fi, q, Aq) and calculate the LOFF state free energy as a function of 6fi and q, for 
given coupling G, average chemical potential fl, and cutoff A. We vary q to minimize 
the LOFF free energy, and compare it with that for the standard BCS pairing ( |4.10D 
to see which is favored. In this way we can map out the phase diagram for the three 
phases of pairing between the two species of quark: BCS, LOFF, and unpaired. 
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Note that the solution to the gap equation, the LOFF gap parameter A^, is not 
a gap in the spectrum of excitations. The quasiparticle dispersion relations ( [4 .71) 
vary with the direction of the momentum, yielding gaps that vary from zero (for 
momenta on the edge of the blocking regions in phase space) up to a maximum of 

We will first discuss the range of in which there exists a LOFF state as a local 
energy minimum. Later we will go on to study the competition between LOFF and 
BCS, and see in what range of the LOFF state is the global minimum. We expect 
the BCS state to be preferred when the mismatch 5^ between the Fermi energies 
of the two species is small. When the mismatch is comparable to the BCS gap 
(5yU ~ Ao) we expect a transition to LOFF, and at larger we expect all pairing 
to cease. These expectations are largely borne out. 

In general we fix A = 1 GeV and /i = 0.4 GeV, and study different coupling 
strengths G which we parameterize by the physical quantity Aq, the BCS gap of 
Eq. (|4.10| ) which increases monotonically with increasing G. When we wish to 
study the dependence on the cutoff, we vary A while at the same time varying 
the coupling G such that Aq is kept fixed. (This is in the same spirit as using a 
renormalization condition on a physical quantity — Aq — to fix the "bare" coupling — 
G.) We expect that the relation between other physical quantities and Aq will be 
reasonably insensitive to variation of the cutoff A. 

We wish to determine 5/^2, the boundary separating the LOFF phase and the 
normal phase. The LOFF to unpaired phase transition is second order, so it occurs 
where the solution A^ to the LOFF gap equation ( |4.y| ) is zero. Setting A^i = in the 
gap equation (|4.1CI| ) yields an analytical expression relating and g, for any given 
G and A. In Figure ^ we show the A^ = curve for three couplings corresponding 
to Ao = 0.1 GeV (strong couphng), Aq = 0.04 GeV and Aq = 0.01 GeV (weak 
coupling). We have only drawn the zero-gap curve in the region where q > 6^. We 
expect this to be the region of interest for LOFF pairing because when q > Sfi the 
two spheres of Figure ^ do in fact intersect. We have verified that, as described in 
some detail in Ref. |T3|, there are regions of Figure ^a with q < 6fi within which 
the LOFF gap equation (|49|) has (one or even two) nonzero solutions, but these 
solutions all correspond to phases whose free energy is either greater than that of 
the normal phase or greater than that of the BCS phase or both. Figure |^ shows 
that for a given coupling strength, parameterized by Aq, there is a maximum 
for which the LOFF state exists: we call it Sfi2- For 6fi > 6fi2, the mismatch of 
chemical potentials is too great for the LOFF phase to exist. 

We see from Figure |^a that as the coupling gets weaker, 61^2/^0 gets gradually 
larger. (Of course, 6fi2 itself gets smaller: the quantities plotted are Sfi/Ao and 
g/Ao.) Note that in the Aq — >■ limit, the zero gap curve is essentially that shown 
in the figure for Aq = 0.01 GeV, in agreement with the curve obtained at weak 
coupling by Fulde and Ferrell |]T3. The fact that this curve ceases to move in the 



Ao ^ limit means that 5/i2 — > while 8^2/^0 const in this limit. 
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Figure 3: (a) The zero-gap curves for the LOFF state. To the right of a sohd curve, 
there is no solution to the LOFF gap equation, to the left of the curve there is a 
solution, and on the curve the gap parameter is zero. The three curves are (from 
strongest to weakest coupling): Aq = 0.1,0.04,0.01 GeV. The region q < 6fi is 
complicated to describe |0 , and solutions found in this region never give the lowest 
free energy state at a given 6fi. (b) Here, we choose Aq = 0.04 GeV and focus on 
the region near (5/X2, the maximum value of 5// at which the LOFF state exists. The 
dashed curve shows the value of |q| which minimizes the free energy of the LOFF 
state at a given Sfii, discussed below, is also indicated. 



For 6fi —>■ from below, we see from Figure |^ that there is a solution to the 
LOFF gap equation only at a single value of q. For example, at Aq = 0.04 GeV we 



find q = O.SSOAq = 1.183 5/i2 at 5/i2 = 0.744Ao. (In agreement with Refs. pi, p 



in the weak coupling limit we find q = 0.906Ao = 1.20 5/i2 at 6fi2 = 0.754Ao.) For 
any value of 5/i < 5/^2, solutions to the LOFF gap equation exist for a range of |q|. 
We must now find the value of |q| for which the free energy of the LOFF state is 
minimized. We obtain the free energy of the LOFF state at a point in Figure ^ by 
first solving the gap equation ( [4. 91) numerically to obtain A^, and then using ( [4.4| ) 
and ( [4.6|) to evaluate {Fq + Hj) given in ( [4.2|) and (|4.3|) . For each value of 6fi < 6^2 
we can now determine which choice of q yields the lowest free energy. The resulting 
"best-g curve" curve is shown in Figure for Aq = 0.04 GeV.0 

Finally, for each point on the best-g curve we ask whether the LOFF free energy 

^As a check on our determination of the best we have confirmed that the total momentum of 
the LOFF state with the best q is zero, as must be the case for the ground state of the system at a 
given (5/i (by a theorem attributed to Bloch j5^). This is a powerful check, because it requires the 
net momentum of the unpaired quarks in the blocking regions (which is in the negative z direction; 
see Figure ||) to be cancelled by the net momentum carried by the LOFF condensates. When, in 
future work, our ansatz is extended to describe a LOFF crystal rather than a single plane wave, 
this check will no longer be powerful. Once we go from F ~ exp(2iq • r) to F ^ cos(2q • r) or to 
a more involved crystalline pattern, the total momentum of the condensates and of the unpaired 
quarks will each be zero. 
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Figure 4: LOFF and BCS gaps and free energies as a function of 5ji, with coupling 
chosen so that Aq = 40 MeV and with Ji = 0.4 GeV, A = 1 GeV. Free energies 
are measured relative to the normal state. At each we have varied q to find the 
best LOFF state. The vertical dashed line marks 6fi = 6fii, the value of above 
which the LOFF state has lower free energy than BCS. The expanded inset (wherein 
s = 10^^ GeV^) focuses on the region < 5fx < where the LOFF state has 
the lowest free energy. 



at that 5fi and (best) q is more or less than the free energy of the BCS state 
at the same In this way, we find 6ni at which a first order phase transition 
between the LOFF and BCS states occurs. In Figure ^ we show the competition 
between the BCS and LOFF states as a function of the Fermi surface mismatch 
6fi, for a fixed coupling corresponding to Aq = 40 MeV. The LOFF state exists for 
6fi < 6fi2 = 0.744Ao. At each 6fi < 6fi2, we plot the gap parameter and free energy 
characterizing the LOFF state with the best q for that Sfi. Although the BCS gap 
Aq is larger than the LOFF gap A^, as Sfi increases we see from Eq. ( p. 3D that the 
BCS state pays a steadily increasing free-energetic price for maintaining = Pf, 
whereas the LOFF state pays no such price. We now see that the LOFF state has 
lower free energy than the BCS state for 6fi > in this case = 0.7104Ao. At 
6fi = 5 III, the gap parameter is A^ = 0.0078 GeV = O.IQSAq. (Had we calculated 
5^1 by comparing the BCS free energy with that of the unpaired state instead of 
with that of the LOFF state, we would have obtained = 0.711Ao. As the inset 
to Figure ^ confirms, the BCS free energy varies so rapidly that this makes an almost 
imperceptible difference. In later figures, we therefore obtain via the simpler 
route of comparing BCS vs. normal.) At the coupling corresponding to Aq = 40 
MeV, we have found that the LOFF state is favored over both the BCS state and 
the normal state in a "LOFF window" 0.710 < 5/i/Ao < 0.744. 

With solutions to the gap equation in hand, we can obtain the J = condensate 
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Figure 5: The two LOFF condensates (J = 0) and F^ (J = 1) for the same 
choice of parameters as in Figure We focus on the region < 6fi < 6fi2. For 
reference, in the BCS phase F^ = Aq/G = 0.00583 GeV^ and F^ = 0. 



Ta = GAa and the J = 1 condensate F^ given in Eq. ( 4.11|) . In Figure ^, we show 



both condensates within the LOFF window 5fii < 6fj, < 6fj,2- We see first of all that 
Fb 7^ 0, as advertised. For the choice of parameters in Figs. § and || we find F^/F^ 
essentially constant over the whole LOFF window, varying from 0.121 at 6fii to 
0.133 at 6fi2- Increasing Aq tends to increase Tb/Ta, as does decreasing A. Second 
of all, we see that the phase transition at 5/x = Sn2, between the LOFF and normal 
phases, is second order in the mean- field approximation we employ throughout. 

Near the second-order critical point 6fi2, we can describe the phase transition 
with a Ginzburg-Landau effective potential. The order parameter for the LOFF-to- 
normal phase transition is 

$(r) = -^(6,,6,^3#°(r)C75^^-^(r)) (5.1) 

so that in the normal phase $(r) = 0, while in the LOFF phase $(r) = F^e*^*''". 
Expressing the order parameter in terms of its Fourier modes $(k), we write the 
LOFF free energy (relative to the normal state) as 

F{{m}) = E Ic2{e)\m\' + c.ik'mk)]' + om^)) . (5.2) 

k 

For Sfi > 6fi2, all of the C2{k'^) are positive and the normal state is stable. Just 
below the critical point, all of the modes $(k) are stable except those on the sphere 



21 



|k| = 2^2, where q2 is the value of |q| at 5^2 (so that q2 — 1.2Sfi2 — O.QAq at weak 
coupling). In general, therefore, many modes on this sphere can become nonzero, 
giving a condensate with a complex crystal structure. We consider the simplest 
case of a plane wave condensate where only the one mode $(k = 2q2) = Ta is 
nonvanishing. Dropping all other modes, we have 

F{Ta) = a{6fi - 6fi2){TAy + KTa)^ (5.3) 

where a and b are positive constants. Finding the minimum-energy solution for 
6fi < 6fi2, we obtain simple power-law relations for the condensate and the free 
energy: 

r^(5/i) = Kr{6fi2 - 5yu)'/', F{6fi) = -Kf{6^12 - 5/i)'. (5.4) 

These expressions agree well with the numerical results shown in Figs. ^ and ^. The 
Ginzburg-Landau method does not specify the proportionality factors Ky and Kp, 
but analytical expressions for these coefficients can be obtained in the weak coupling 
limit by explicitly solving the gap equation [^, yielding 



GaKt = 2v^V(?2/5/i2)2 - 1 ^ 1.15x/A^ 

Kf = {Aji^/Ti^){{q2/5ii2? - 1) ^ ^ ' ' 



Notice that because {6112 — ^l^i)l^l^2 is quite small, the power-law relations ( WM are 
a good model of the system throughout the entire LOFF interval < 6fi < 8112 
where the LOFF phase is favored over the BCS phase. The Ginzburg-Landau 
expression ( p.3| ) gives the free energy of the LOFF phase near 8^2, but it cannot 
be used to determine the location of the first-order phase transition where the 
LOFF window terminates (locating the first-order point requires a comparison of 
LOFF and BCS free energies). 

It is interesting to explore how the width of the LOFF window depends on the 
strength of the coupling, and to confirm that it is insensitive to the cutoff. We do 
this in Figure H, where we plot 6^12/^0 (solid lines) and 5/xi/Ao (dashed fines). The 
LOFF state is favored for 6fii/AQ < 6fi/Ao < 6^12/^.0, i.e. between the sofid and 
dashed curves in Figure ^ In the weak coupling limit, the LOFF window tends 
to 0.707 < (5/i/Ao < 0.754 and at tends to 0.23Ao, as in Refs. [|T|, |T§. 
Note that if one takes the weak-coupling limit Aq ^ at fixed neither BCS 
nor LOFF pairing survives because S/i/Aq — * 00. However, for any arbitrarily small 
but nonzero coupling, the LOFF phase is favored within a range of 6fi. Figure ^ 
thus demonstrates that in an analysis of the LOFF state in the weak-coupling limit, 
it is convenient to keep S^i/Aq fixed while taking Aq — > 0. We see from Figure | 
that strong coupling helps the BCS state more than it helps the LOFF state. When 
the coupling gets strong enough, there is no longer any window of Fermi surface 
mismatch 6fi in which the LOFF state occurs: the BCS state is always preferred. 

The different lines of each type in Figure |^ are for different cutoffs and show 
that there is in fact little sensitivity to the cutoff. The A dependence of 6^1/ Aq and 
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Figure 6: The interval of 5/i within which the LOFF state occurs, as a function of the 
couphng (parameterized as usual by the BCS gap Aq). Below the solid line, there 
is a LOFF state. Below the dashed line, the BCS state is favored. The different 
lines of each type correspond to different cutoffs A = 0.8 GeV to 1.6 GeV. 5/ii/Ao 
and (5/i2/Ao show little cutoff-dependence, and the cutoff-dependence disappears 
completely as Aq, 5/i ^ 0. 

(5/i2/Ao is mild for all values of Aq which are of interest, and is weakest for Aq 0. 
This is because in that limit pairing can only occur very close to the unblocked 
ribbon of the ellipsoid of Fig. 2b, along which the integrand in the gap equation is 
singular and pairing is allowed. Thus most of the pairing region V, and in particular 
the region near A, become irrelevant in this limit. 

The one physical quantity which we have explored which does turn out to depend 
qualitatively on A is the ratio Tb/^a- Those quarks with momenta as large as A 
which pair have momenta which are almost antiparallel, and so contribute much 
less to than to F^. For this reason, the ratio Tb/^a is sensitive to the number 
of Cooper pairs formed at very large p, and hence to the choice of A. As discussed 
above, pairing far from the favored ribbon in phase space becomes irrelevant for 
Aq 0, and indeed in this limit we find that the A dependence of Tb/Ta decreases. 
However, for Aq = 40 MeV we find that changing A from 1.2 GeV to 0.8 GeV 
increases Tb/Ta by more than 50%. 

We chose to show results for Aq = 40 MeV in Figure ^ because with this choice, 
the LOFF window occurs at values of 6fi comparable to that in the illustrative 
example Q: 6fi = - = 27 MeV. Of course, neither 5/i nor the value of 
Aq are accurately known for the quark matter which may exist within a compact 
star. Still, it seems possible that their ratio could be appropriate for the quark 
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matter to be in the LOFF phase. If there is a range of radii within a compact star 
in which quark matter occurs with 6fii < 6^ < 6fi2, this quark matter will be a 
crystalline color superconductor. 

In Figure ^ the LOFF gap parameter is 7.8 MeV at 6fi = 6fii. It remains 
larger than typical neutron star temperatures ~ 1 keV until very close to 6fi = 
(5/i2. Similarly, the LOFF free energy, which is 4.8 x 10"*^ GeV^ = 4.8 x (10 MeV)^ 
at 6fi = 6fii, is much larger than T^^ throughout the LOFF window except very 



close to 5/i = 6fi2- Furthermore, we shall see in Section 7.3 that the free energy of 



the LOFF state is of the right order to lead to interesting glitch phenomena. 
6 More general Hamiltonian and ansatz 

In Section ^ we introduced the four-fermion interaction Hamiltonian Hj of Eq. ( [4.1D 
with independent couplings Ge and Gm for the interactions which model the 
exchange of electric and magnetic gluons. It proves convenient to use the linear 
combinations 

Ga = j{Ge + ^Gm) ^.^ 
Gb = j{Ge — Gm) , 

of the coupling constants in terms of which the expectation value of Hj in the LOFF 
state ( |3.2| ) becomes 

{Hj) = -\GaV (ir^p + |r«p) - \GsV (irip + |r|p) . (6.2) 

Thus, a positive coupling Ga describes an attractive interaction which induces a 
J = condensate F^. As we have seen, in the LOFF state this is necessarily 
accompanied by a J = 1 condensate F^. In our analysis to this point, we have set 
Ga = G > and = 0. We now discuss the general case, in which 7^ 0. 

Before beginning, let us consider how to choose Gb/Ga in order for our model 
Hamiltonian to be a reasonable toy model for QCD at nonzero baryon density. At 
zero density, of course, Lorentz invariance requires = 0. At high densities, on 
the other hand, electric gluons are screened while static magnetic gluons are not. 



(Magnetic gluons with nonzero frequency are damped.) We now know |[T^ that at 
asymptotically high densities it is in fact the exchange of magnetic gluons which 
dominates the pairing interaction. This suggests the choice G^; = 0, corresponding 
to Gb/Ga = —1/3. At the accessible densities of interest to us, it is presumably not 
appropriate to neglect Ge completely. Note also that the four-fermion interaction 
induced by instantons in QCD only yields interactions in flavor-antisymmetric 
channels. It results in an attractive interaction in the J = channel and no 
interaction in the J = 1 channel. Thus, although the instanton interaction cannot 
be written in the form ( [4.1|) , for our purposes it can be thought of as adding a 



contribution to Ga-, but none to Gb. Hence our model is likely to best represent 
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high density QCD for a ratio of couphngs lying somewhere in the range 

-i<^<0. (6.3) 

We plot our results over a wider range of couplings below. 

Once Gb ^ and there is an interaction in the J = 1 channel, we expect, in 
addition to the J = 1 condensate F^, a J = 1 gap parameter A^. The quasiparticle 
dispersion relations are then determined by A a and A^, which are defined as 

A^ = GaTa 4^ 
Ab — Gb^b ■ 

Following through the variational calculation as in Section || leads to the coupled 
gap equations: 

^ _ 2(7^ ^ 25'a(Aa5'a + AbSb) 

^ ~ ~^^^/{h + p\ + \q-p\~2fiy + 4{AASA + A^ 
^ _ '2Gb sr-^ 2Sb{AaSa + AbSb) 

~ ^ + pI + lq - pI - 2/^)' + 4(Aa5a + AbSbY ^^-^^ 

Sa = sin(i/3yi(p)) 
Sb = sin(i/5B(p)) 

with /3a(p) = a«(p) + «d(p), Pb{p) = ^^(p) — ad(p) defined in terms of the 
angles in Figure |I|. The pairing region V is still defined by ( |4.8| ) but with new 
quasiparticle dispersion relations obtained from Eqs. ( [4.7|) with A\Sa replaced by 
(AaSa + AbSb?. 

For Gb = 0, the coupled equations ( |6.5|) reduce to Eqs. ( [4.9| ) and ([4.11| ). Note 



that if, instead, Gb > and Ga = 0, we find an attractive interaction in the J = 1 
channel in Eq. ( |6.1|) and no interaction in the J = channel. Analysis of Eqs. ( |6.5|) 
in this case yields a nonzero value of A^, while A a = even though F^ 7^ 0. The 
geometry of the LOFF pairs requires Ta ^ when F^ 7^ 0. 

Rather than describing how every Figure in Section |^ changes when Gb ^ 0, we 
choose to focus on the question of how the interval of Sfi within which the LOFF 
state occurs (the LOFF window) changes as a function of Gb/Ga- To further 
simplify the presentation, we specialize to the weak-coupling limit in which Aq —>■ 0. 
This means that, as in Figure ^, the LOFF window is independent of the cutoff A. 

We show the dependence of the LOFF window on Gb/Ga in Figure]^. The lower 
boundary 6fi = 6fi^ is, as in Section |^, the same (up to a very small correction) as 
the 6fi at which the BCS and normal states have equal free energies. We find the 
upper boundary 6fi = 6fi2 by first dividing Eqs. (|6.5| ) by Aa and then looking for a 
value of 6fi at which A^ and A^ ^ but Aa/Ab remains nonzero. As before, 
this defines a zero-gap curve, and Sfi2 is the maximum value of 6^ reached by this 
curve. 
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Figure 7: The interval of 5^ in which the LOFF state is favored at weak couphng, 
as a function of the ratio of couphngs Gb/Ga- Below the solid line, there is a LOFF 
state. Below the dashed line, the ordinary BCS state is favored. Gb = corresponds 
to the Lorentz-invariant interaction with Ge = Gm- QCD at high density is likely 
best described by a coupling in the range — | < Gb/Ga < 0. 




We find that the lower boundary d^ii is completely unaffected by the value of 
Gbi since the BCS state is purely J = 0. So in the weak-coupling limit we obtain 
the result of Section]^, 5yU^/Ao = 0.707, independent of Gb/Ga- In contrast, ^/ij, 
the upper boundary of the LOFF window, increases with increasing Gb- This is 
understandable: the LOFF state always produces a J = 1 condensate, so we expect 
it to be fortified by > and penalized by < 0. There is no analogue of 
this behavior in an electron superconductor [^, where there can be no J = 1 
condensate. Our J = 1 condensate affects the gap equation and free energy only if 
7^ 0; for this reason, our weak coupling results are in agreement with those of 
LOFF |]T2|, |T3| only if G^ = 0, as in Section ||. The effect of a coupling G^ in the 
physically interesting range (|6.3| ) is to reduce the LOFF window, but only slightly. 

In both this Section and the previous one, we have calculated by examining 
the competition between LOFF pairing and no pairing. Should we instead have 
considered the competition between LOFF pairing and the formation of {uu) and 
{dd) condensates, each at their respective Fermi surface, each with q = 0? Assuming 
as usual that the color antisymmetric channel is the most attractive one, flavor 
symmetric pairing requires spin-symmetric pairing, i.e. J = 1 [Q. Within the 
model and ansatz that we have considered, the question is easily answered. If we 
choose Ga > and Gb = 0, as in Sections ^ to ^, there is no interaction in the 
spin-symmetric, flavor-symmetric, color-antisymmetric channel. If we strengthen 
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magnetic gluon exchange relative to electric gluon exchange by choosing < 0, 
the interaction in this channel is repulsive. We have confirmed this by evaluating the 
expectation value of Hi in a state with spatially uniform J = 1 pairing and (uCa^^u) 
condensate (obtained by using two u creation operators in the ansatz (|3.2|) , setting 
= as appropriate for (uu) pairing, setting q = and removing the blocking 
regions). We find that Ga gives no interaction in this channel and Gb < is 
repulsive. Thus, for the same reason that the LOFF window shrinks for < 0, 
there can be no (uu) or (dd) pairing. However, the scenario is apparently different 
at asymptotically high density: it has been shown by Schafer that long-range 



single-gluon exchange does in fact induce pairing in this J = 1 channel. (The long- 
range interaction emphasizes near-coUinear scattering which is attractive for both 
electric and magnetic gluons.) For either a pointlike interaction with > or a 
long-range interaction dominated by near-collinear scattering, we therefore expect 
competition between LOFF pairing and (uu) and (dd) pairing, since the latter would 
then be favored for 6fi > 6^2- 

Our ansatz only contains LL and RR pairing. We leave a complete analysis of 
the generalization to LR pairing to future work. We have, however, constructed the 
ansatz for spatially uniform LR pairing with = and q = 0. We find that the 
interaction in this J = 1 channel is attractive if G^ + Ga/ > and is independent 
of the linear combination of couplings Ge — ^Gm- The J = 1 channel with LR 
pairing yields a {uC'^^u) condensate instead of the J = 1 condensate (uCa^'^u) 



obtained for the case of LL and RR pairing. In agreement with Ref. we find 
that magnetic gluon exchange, with Ge = and Gm > 0, is attractive in the 
(uCYu) channel. Note that in the nonrelativistic limit {uC'^'^u) and {uCa^^u) are 
equivalent J = 1 condensates. In the relativistic setting relevant in quark matter, 
we find that pointlike interactions in these two channels have opposite sign. 

We have set up the gap equation describing a spatially uniform {uCYu) 
condensate and solved it for G^; = Gm = G, Aq = 40 MeV, /i„ = 0.4 GeV, 
and A = 1 GeV. We find a gap of 8 keV and a free energy which is about five orders 
of magnitude smaller than that of the LOFF phase. (If we choose G^; = and 
Gm > 0, the interaction is still attractive but the gap is even smaller.) Therefore, 
even though for 6fi > 5^2 we expect LR pairing and consequent {uCYu) and 
(dGYd) condensates, the resulting condensation energy is so small that it is a good 
approximation to neglect these condensates in the evaluation of Sfi2, as we have 
done. We leave for future work a complete analysis of the competition between the 
LOFF phase (with an ansatz extended to allow Li?-LOFF pairing) and the spatially 
uniform (mG7*m) condensate. 



7 Conclusions, future work and astrophysical implications 
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7.1 Conclusions 



We have studied the formation of a rotational-symmetry-breaking LOFF state 
involving pairing between two flavors of quark whose chemical potentials differ by 
2Sii. This state is characterized by a gap parameter and a diquark condensate, but 
not by an energy gap in the dispersion relation. In the LOFF state, each Cooper 
pair carries momentum 2q with |q| ^ 1.25/i. The condensate and gap parameter 
vary in space with wavelength vr/|q|. 

We focused primarily on an NJL-type four-fermion interaction with the quantum 
numbers of single gluon exchange. In the limit of weak coupling (BCS gap Aq ^ yu) 
the LOFF state is favored for values of 5/i which satisfy 6fii < 6fi < 5/i2, where 
Sfi^/Ao = 0.707 and 6fi2/^o = 0.754. The LOFF gap parameter decreases from 
0.23Ao at 6fi = 6fii to zero at 6fi = 5^2- These are the same results found by LOFF 
in their original analysis. Except for very close to 5/i2, the critical temperature 
above which the LOFF state melts will be much higher than typical neutron star 
temperatures. At stronger coupling the LOFF gap parameter decreases relative to 
Aq and the window of (5/i/Ao within which the LOFF state is favored shrinks. The 
window grows if the interaction is changed to weight electric gluon exchange more 
heavily than magnetic gluon exchange. 

Because it violates rotational invariance by involving Cooper pairs whose 
momenta are not antiparallel, the quark matter LOFF state necessarily features 
nonzero condensates in both the J = and J = 1 channels. Both condensates are 
present even if there is no interaction in the J = 1 channel. In this case, however, 
the J = 1 condensate does not affect the quasiparticle dispersion relations; that is, 
the J = 1 gap parameter vanishes. If there is an attraction in the J = 1 channel 
(as, for example, if the strength of the electric gluon interaction is increased) the 
size of the LOFF window increases. 

The quark matter which may be present within a compact star will be in the 
crystalline color superconductor (LOFF) state if Sfi/Ao is in the requisite range. For 
6fi as in the illustrative example ( |2.1| ), this occurs if the gap Aq which characterizes 
the uniform color superconductor present at smaller values of 6fi is about 40 MeV. 
This is in the middle of the range of present estimates. Both 6fi and Aq vary as a 
function of density and hence as a function of radius in a compact star. Although it 
is too early to make quantitative predictions, the numbers are such that crystalline 
color superconducting quark matter may very well occur in a range of radii within 
a compact star. It is therefore worthwhile to consider the consequences. 

7.2 Future work 

The prospect of spontaneous violation of translational and rotational symmetry 
in dense quark matter is very exciting. In the final pages of this paper we will 
begin to explore one particularly interesting consequence: glitch behavior in quark 
matter within compact stars. First, however, we list a number of direct extensions 
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of our work, several of which are prerequisites to a quantitative exploration of 
the astrophysical consequences of crystalline color superconducting quark matter 
in compact stars. 

1. We have restricted ourselves to two flavors of quark, and varied freely. It 
is crucial to look at more realistic examples, imposing charge neutrality and 
weak equilibrium, and including the strange quark. We expect a LOFF phase 
wherever {us), (ds) or (ud) pairs approach their unpairing transitions, but 
this must be verified quantitatively. Further generalizations would include 
bare quark masses and spontaneous generation of constituent quark masses 
by chiral condensation. 

2. It would be valuable to complement our NJL-model study with a controlled 
calculation using one gluon exchange in the asymptotically high density limit. 
There are two reasons why this is worthwhile. First, it will allow a controlled 
analysis without model assumptions, albeit one of quantitative value only at 
extremely high densities. In particular, this would allow a better estimation 
of the relative magnitude of the J = 1 and J = condensates, which was the 
one feature which we found to depend strongly on the choice of cutoff in our 
model. Second, quark-quark scattering by the exchange of a gluon at weak 
coupling is dominated by small-angle scattering, whereas in an NJL model 
of the type we have used this is not the case. This can actually affect the 
sign of the interaction in the J = 1 channel and perhaps thereby increase 
the range of the LOFF window, as we pointed out in the previous section. 
Moreover, it is known that the LOFF window is much wider in one dimension 



than in three and since the three-dimensional physics at asymptotically 
high densities can be treated as a sum of one-dimensional theories |]2T| , |26| , 
we have another reason to suspect that the LOFF window may be wider at 
asymptotically high densities than our present analysis would suggest. 

3. As we have discussed at length in Section H, it would be of interest to extend 
our treatment to include pairing between quarks of the same flavor and pairing 
between quarks of opposite chirality. 

4. Perhaps the most crucial unresolved issue is the question of what crystal 
structure the LOFF phase chooses. Larkin and Ovchinnikov concluded that 
the condensate varies in space like cos(2q-r), forming a one-dimensional 
standing wave with nodal planes spaced every 7r/(2|q|). The competition 
between this planar structure and one with, say, a cubic or body-centered-cubic 
crystal structure is subtle. In two dimensions, the answer depends sensitively 



on the temperature [^]; in three dimensions, it is apparently still unresolved 



even in the original LOFF context [H2|. In the QCD context, with the added 



complication of a J = 1 condensate, it will be quite interesting to determine 
what pattern is favored. 
5. Finally, it would be very interesting to investigate the astrophysical 
consequences of the LOFF phase. Since it occurs in a range of 6fi, one would 
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expect that quark matter stars could contain a layer of crystalline LOFF 
condensate. In the next subsection, we take some preliminary steps in this 
investigation. 



7.3 Looking ahead to astrophysical consequences 

Many pulsars have been observed to glitch. Glitches are sudden jumps in rotation 
frequency f2 which may be as large as Ail/n ~ 10"^ but may also be several orders 
of magnitude smaller. The frequency of observed glitches is statistically consistent 
with the hypothesis that all radio pulsars experience glitches |5^. Glitches are 
thought to originate from interactions between the rigid crust, somewhat more 
than a kilometer thick in a typical neutron star, and rotational vortices in the 
neutron superfiuid. The inner kilometer of the crust consists of a rigid lattice 
of nuclei immersed in a neutron superfiuid |6y]. Because the pulsar is spinning, 



the neutron superfiuid (both within the inner crust and deeper inside the star) is 
threaded with a regular array of rotational vortices. As the pulsar's spin gradually 
slows due to emission of electromagnetic radiation, these vortices must gradually 
move outwards since the rotation frequency of a superfiuid is proportional to the 
density of vortices. Deep within the star, the vortices are free to move outwards. 
In the crust, however, the vortices are pinned by their interaction with the nuclear 
lattice. What happens next varies from model to model. Perhaps the vortices exert 
sufficient force on the crust to tear it apart, resulting in a sudden breaking and 
rearrangement of the crust and a change in the moment of inertia [^. Perhaps a 



large cluster of vortices within the inner crust builds up enough outward pressure to 
overcome the pinning force, suddenly becomes unpinned, and moves macroscopically 



outward |£2| |6^, ^ |6§, p^, This sudden decrease in the angular momentum 
of the superfiuid within the crust results in a sudden increase in angular momentum 
of the rigid crust itself, and hence a glitch. Perhaps, due to interactions between 
neutron vortices and proton fiux tubes, the neutron vortices pile up just inside the 
inner crust before suddenly coming unpinned [^. Although the models differ in 
important respects, all agree that the fundamental requirements are the presence 
of rotational vortices in a superfiuid and the presence of a rigid structure which 
impedes the motion of vortices and which encompasses enough of the volume of the 
pulsar to contribute significantly to the total moment of inertia.0 

Although it is premature to draw quantitative conclusions, it is interesting to 
speculate that some glitches may originate deep within a pulsar which features a 
quark matter core, in a region of that core in which the color superconducting quark 
matter is in a LOFF crystalline color superconductor phase. The first prerequisite 
for a quantitative answer to whether this may occur is to repeat our analysis in 



^The first model of glitches which was proposed |7C[] relies on the cracking and settling of the 
neutron star crust ( "starquakes" ) as the neutron star spins down. This model does not require the 
presence of rotational vortices. However, this model fails to explain the magnitude and frequency 
of glitches in the Vela pulsar [E^ ^sl . 
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the more general context of three-flavor quark matter with a nonzero strange quark 
mass Ms, to estimate over what range of densities LOFF phases may arise, as either 
{ud), {us) or {ds) condensates approach their unpairing transitions. Comparison 
to existing models which describe how p^, p'^ and pp vary within a quark matter 
core in a neutron star would then permit an estimate of how much the LOFF 



region contributes to the moment of inertia of the pulsar. Furthermore, a three 
flavor analysis is required to determine whether the LOFF phase is a superfluid. If 
the only pairing is between u and d quarks, this 2SC phase is not a superfluid 0, |[, 
whereas if all three quarks pair in some way, a superfluid is obtained 0, |] .[] 

Henceforth, we suppose that the LOFF phase is a superfluid, which means that 
if it occurs within a pulsar it will be threaded by an array of rotational vortices. It is 
reasonable to expect that these vortices will be pinned in a LOFF crystal, in which 
the diquark condensate varies periodically in space. Indeed, one of the suggestions 
for how to look for a LOFF phase in terrestrial electron superconductors relies on 
the fact that the pinning of magnetic flux tubes (which, like the rotational vortices 
of interest to us, have normal cores) is expected to be much stronger in a LOFF 
phase than in a uniform BCS superconductor |^ . 



A real calculation of the pinning force experienced by a vortex in a crystalline 
color superconductor must await the determination of the crystal structure of the 
LOFF phase. We can, however, attempt an order of magnitude estimate along 



the same lines as that done by Anderson and Itoh for neutron vortices in the 
inner crust of a neutron star. In that context, this estimate has since been made 
quantitative ||7^, BB, With parameters chosen as in Figure ^ we find that 



at 5/i = the LOFF phase is favored over the normal state by a free energy 
-Floff ~ 5 X (10 MeV)'^ and the spacing between nodes in the LOFF crystal is 
6 = 7r/(2|q|)~9 fm. The thickness of a rotational vortex is given by the correlation 
length ^ ~ 1/A^ ~ 25 fm. All these numbers are quite uncertain, but we will use 
them for the present. In the context of crustal neutron superfluid vortices, there 
are three distinct length scales: the vortex thickness ^, the lattice spacing between 
nuclei 6, and i?, the radius of the individual nuclei. (The condensate vanishes within 
regions of size R separated by spacing h.) In the LOFF phase, the latter two length 
scales are comparable: since the condensate varies like cos(7rr/6) it is as if i? ~ h. 
The fact that these length scales are similar in the LOFF phase will complicate 
a quantitative calculation of the pinning energy; it makes our order of magnitude 
estimation easier, however. The pinning energy is the difference between the energy 
of a section of vortex of length h which is centered on a node of the LOFF crystal 
vs. one which is centered on a maximum of the LOFF crystal. It is of order 

Ep ~ Floff ~ 4 MeV . (7.1) 



^As an aside, note that the crystalhne chiral condensate (due to particle-hole pairing which 
may form at sufficiently strong coupling or at very large Nc) is not a superfluid. 
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The resulting pinning force per unit length of vortex is of order 



f 4MeV 

A complete calculation will be challenging because b < C,, and is likely to yield an fp 
which is somewhat less than that we have obtained by dimensional analysis |^6|, |68[ . 
Note that our estimate of fp is quite uncertain both because it is only based on 
dimensional analysis and because the values of A a, b and -Floff are uncertain. (We 
know the values of all the ratios A^/Aq, (5/i/Ao, Q'/Aq and consequently 6Ao quite 
accurately in the LOFF phase. It is of course the value of the BCS gap Aq which is 
uncertain.) It is therefore premature to compare our crude result to the results of 
serious calculations of the pinning of crustal neutron vortices as in Refs. ^ |BS 



It is nevertheless remarkable that they prove to be similar: the pinning energy of 
neutron vortices in the inner crust is 

^ 1 - 3 MeV (7.3) 



and the pinning force per unit length is ||66 



•'^ b^ (25 - 50 fm) (4-20 fm) ' ^ ' ^ 

where the form of this expression is appropriate because < b. Perhaps, therefore, 
glitches occurring in a region of crystalline color superconducting quark matter may 
yield similar phenomenology to those occurring in the inner crust. 

The reader may be concerned that a glitch deep within the quark matter 
core of a neutron star may not be observable: the vortices within the crystalline 
color superconductor region suddenly unpin and leap outward; this loss of angular 
momentum is compensated by a gain in angular momentum of the layer outside the 
LOFF region; how quickly, then, does this increase in angular momentum manifest 
itself at the surface of the star as a glitch? If the LOFF layer is the outer layer 
of the quark matter core — not unreasonable since the chemical potential differences 
will be larger here than deeper inside the quark matter — there is no problem. The 
LOFF glitch speeds up the nucleon superfluid outside the quark matter core, and 
the rotation of this superfluid is coupled to the rotation of the outer crust on very 
short time scales . This rapid coupling, due to electron scattering off vortices and 



the fact that the electron fluid penetrates throughout the star, is usually invoked to 
explain that the core nucleon superfluid speeds up quickly after a crustal glitch: the 



only long relaxation time is that of the vortices within the inner crust |T^. Here, we 
invoke it to explain that the outer crust speeds up rapidly after a LOFF glitch has 
accelerated the quark matter at the base of the nucleon superfluid. After a glitch 
in the LOFF region, the only long relaxation times are those of the vortices in the 
LOFF region and in the inner crust. 
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A quantitative theory of glitches originating within quark matter in a LOFF 
phase must await the further microscopic calculations sketched in Section |7.2| . 
In particular, an understanding of points 1 and 4 of Section |7.2| is a mandatory 
prerequisite. However, our rough estimate of the pinning force on rotational vortices 
in a LOFF region suggests that this force may be comparable in magnitude to that 
on vortices in the inner crust of a conventional neutron star, which yields glitches in 
accord with those observed in pulsars. This is surely strong motivation for further 
investigation. 

Perhaps the most interesting consequence of these speculations arises in the 
context of compact stars made entirely of strange quark matter. The work of 
Witten 1 74] and Farhi and Jaffe [^] raised the possibility that strange quark matter 
may be energetically stable relative to nuclear matter even at zero pressure. If 
this is the case it raises the question whether observed compact stars — pulsars, for 
example — are strange quark stars []7B|, rather than neutron stars. A conventional 
neutron star may feature a core made of strange quark matter, as we have been 
discussing above. Strange quark stars, on the other hand, are made (almost) entirely 
of quark matter with either no hadronic matter content at all or with a thin crust, of 
order one hundred meters thick, which contains no neutron superfluid |l77| , [78| . The 
nuclei in this thin crust are supported above the quark matter by electrostatic forces; 
these forces cannot support a neutron fluid. Because of the absence of superfluid 
neutrons, and because of the thinness of the crust, no successful models of glitches 
in the crust of a strange quark star have been proposed. Since pulsars are observed 
to glitch, the apparent lack of a glitch mechanism for strange quark stars has been 



the strongest argument that pulsars cannot be strange quark stars |^ Q |8T|] . This 
conclusion must now be revisited. 

Madsen's conclusion that a strange quark star is prone to r-mode instability 
due to the absence of damping must also be revisited, since the relevant fluid 
oscillations may be damped within or at the boundary of a region of crystalline 
color superconductor. 

The quark matter in a strange quark star, should one exist, would be a color 
superconductor. Depending on the mass of the star, the quark number densities 
increase by a factor of about two to ten in going from the surface to the center |77 . 



This means that the chemical potential differences among the three quarks will 
vary also, and there could be a range of radii within which the quark matter is in 
a crystalline color superconductor phase. This raises the possibility of glitches in 
strange quark stars. Because the variation in density with radius is gradual, if a 
shell of LOFF quark matter exists it need not be particularly thin. And, we have 
seen, the pinning forces may be comparable in magnitude to those in the inner 
crust of a conventional neutron star. It has recently been suggested (for reasons 
unrelated to our considerations) that certain accreting compact stars may be strange 



quark stars |S2], although the evidence is far from unambiguous [Q. In contrast 



it has been thought that, because they glitch, conventional radio pulsars cannot be 
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strange quark stars. Our work questions this assertion by raising the possibihty that 
ghtches may originate within a layer of quark matter which is in a crystalhne color 
superconducting state. 
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